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Abstract 

The main motivation to investigate hard probes in heavy ion colHsions is to do tomography, 
i.e. to infer medium properties from the in-medium modification of hard processes. Yet while 
the suppression of high Pj hadrons has been measured for some time, sohd tomographic infor- 
mation is slow to emerge. This can be traced back to theoretical uncertainties and ambiguities 
in modelling both medium evolution and parton-medium interaction. Ways to overcome these 
difficulties are to constrain models better and to focus on more differential observables. Correla- 
tions of high Pj hadrons offer non-trivial information beyond what can be deduced from single 
hadron suppression. They reflect not only the hard reaction being modified by the medium, but 
also the back reaction of the medium to the hard probe. Models for hard back-to-back correla- 
tions are now very well constrained by a wealth of data and allow insights into the nature of the 
parton-medium interaction as well as first true tomographic results. Models of full in-medium 
jet evolution are being actively developed, but have yet to make substantial contact with data. 
Progress is slower in the understanding of low Pj correlations, the ridge and the cone, although 
a qualitative understanding of the nature of the physics behind these correlations starts to emerge. 



1. Tomography with hard probes 

In ultrarelativistic heavy-ion (A-A) collisions, a new state of hot and dense QCD matter 
is produced and offers the exciting prospect of studying collectivity and a phase transition in 
a quantum field theory experimentally. The idea to utilize tomography as a tool to study the 
properties of this state of matter is rather simple: The basic principle of tomography is to generate 
a probe with known properties, propagate it through a medium and infer the properties of the 
medium from the changes induced on the probe. In A-A collisions, hard processes occuring 
along with the production of the medium generate suitable tomographic probes in terms of high 
Pj partons. After their initial production with high virtuality, they undergo the quantum evolution 
towards lower virtuality scales while interacting with the surrounding medium, the net result of 
which is the suppression of observed hard processes as compared to the scaled expectation from 
proton-proton (p-p) collisions. 

In order to describe the physics underlying this suppression theoretically, three ingredients 
are needed. Fist, there is the hard process itself, which due to a separation of scales can be treated 
distinct from the medium dynamics and is calculable with good accuracy in perturbative QCD. 
Then, there is the interaction of the evolving parton shower with the medium. This however 
depends on the unknown, non-perturbative properties of the medium (e.g. in terms of relevant 
degrees of freedom) and ultimately needs to be inferred from measurement. Finally, there is the 



Preprint submitted to Nuclear Physics A 



September 11, 2009 



complication that the medium itself evolves at the same timescales set by the propagation of the 
hard probe, thus the evolution of the medium density as a function of space and time needs to 
be known. While there are good constraints from bulk matter observables, the details of this 
evolution are not known with certainty. 

Tomography in the usual sense would be knowing the source and the parton-medium inter- 
action and inferring the density evolution from those. However, since there are two unknowns, 
a different analysis strategy is needed. A promising way to proceed is to focus on all available 
constraints rather than a single high Pj observable as often done. This, in short, is the reason for 
the importance of hard correlation measurements — they offer a handle to observe the medium 
in a different way in terms of in-medium pathlength distribution and type of the high pT parton 
and hence allow ultimately together with other observables to separate the interaction dynamics 
between hard probe and medium from the evolution dynamics of the medium. 

Conceptually, this requires to compute hard processes in a medium model constrained by 
bulk observables such as a hydrodynamical model. In this way, combinations of medium mod- 
els and parton-medium interaction models can be tested against different observables. Thus, 
jet tomography can be a well constrained science if the full wealth of data provided by experi- 
ments is used. An example of such an approach is shown in fP] where three different models for 
the parton-medium interaction (AMY, ASW and HT) are used within the same hydrodynamical 
background to compute hard single hadron observables. It is clearly seen that while the numeri- 
cal extraction of medium parameters leads to values of the transport coefficient q different up to 
a factor two, qualitatively the models are very similar in terms of the Pj dependence of single 
hadron suppression, the extrapolation to different collision centrality or the way they probe ge- 
ometry. To second order however, they predict somewhat different behaviour of the suppression 
with respect to the angle of the hard hadron with the reaction plane, a property which can be 
used for discrimination. In the following, it will be shown that correlation measurements are 
even more powerful to discriminate between models and hence to extract medium properties. 

2. Dihadrons and jets — hard correlations 

Any hard process in QCD creates a back-to-back pair of partons at high virtuality which 
subsequently evolve through a partonic shower with decreasing virtuality into hadron jets. If one 
focuses on single hard hadron observables, one is biased towards events in which a single parton 
carries most of the momentum of the shower. The fragmentation function usually includes effects 
of both perturbative shower evolution and non-perturbative hadronization, but in this particular 
case it chiefly describes the hadronization of the leading parton. Therefore the approximation 
is often made that the medium effect can be cast into the form of energy loss to the medium, 
followed by vacuum fragmentation. The medium-modified fragmentation function for a parton 
with initial energy E and virtuality Q which describes the production of hadrons at momentum 
fraction z from this parton can then be cast into the schematical form 

D„,ed{z, Q) = P(A£, E) ® D,ac{z, Q^Qh)® D,aAz, Qh) (1) 

where P(E, AE) is the probability distribution for energy loss AE in the medium given initial 
energy E and Dyadz, Q Qh) describes the perturbative partonic shower evolution decreasing 
the virtuahty scale from Q down to a hadronic scale Qi, whereas D^.adz, Qh) stands for the non- 
perturbative hadronization process. 
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Figure 1: A schematic picture of tlie surface bias effect for different parametric pathlength dependence of energy loss. 



The energy loss picture is applicable for observables which are sensitive to the leading hadron 
in a shower. This includes single hadron observables and also the hard back-to-back correlation 
with a trigger However, in order to discuss correlations on the same side of the trigger, the 
approximation is not sufficient as it does not treat subleading partons correctly, and in the case 
of fully reconstructed jets where there is no bias towards a single parton carrying most of the 
energy, the energy loss picture is insufficient. Here, the whole parton shower needs to be evolved 
in the medium, whereas hadronization can, due to uncertainty principle arguments, still be as- 
sumed to take place outside the medium. The schematical expression for the medium-modified 
fragmentation function is then 

D„ed(z, Q) = D,„,d(z, Q ^ Qh) ® D,,,(z, Qh). (2) 

2.1. Dihadron correlations and pathlength 

Hard back-to-back hadron correlations probe the leading hadron(s) of each shower and hence 
the energy loss picture can be applied. Calculations have been done in NLO pQCD for a hard 
sphere overlap with Bjorken expansion using the ZOWW energy loss formulation ||^, using LO 
pQCD and the ASW energy loss formulation [3] in 2-d [|4i,|5||] and also as a function of the angle 
with the reaction plane in 3-d hydrodynamical expansions j^. A calculation with elastic energy 
loss in a 3d-hydrodynamics is presented in |7]. Thus, at least some of the systematics is known 
and it can be inferred what properties of the system are reflected in such correlations. 

It emerges from these studies that dihadron correlations are an excellent tool to probe dif- 
ferences in the parametric dependence of energy loss on pathlength. The reason is apparent 
from Fig. [T] If one triggers on a hadron (red) on the near side, the most likely vertex of origin 
will be close to the surface of the medium, and this bias will be stronger for higher powers of 
the pathlength dependence of energy loss. However, this in turn influences the average away 
side pathlength (blue) which becomes larger for higher powers of pathlength dependence and 
moreover is weighted with that higher power Therefore, a distinction between L-pathlength de- 
pendence (in a constant medium) as characteristic for elastic energy loss and dependence is 
very clear in back-to-back correlations and leads to factors two or more difference in observables 
|7]. A conclusion common to the ZOWW |2] and the ASW computation of dihadron correlations 
ILZI] is that the average contribution of elastic energy loss is (9(10%). 

The formalisms differ however in details, such as the geometrical picture of the suppression. 
As shown in Fig. |2] the ZOWW computation shows a strong tangential bias whereas the ASW 
computation does not. Presumably the reason for this difference is that the ZOWW picture uti- 
lizes an averaged energy loss given a path through the medium whereas the ASW picture allows 
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Figure 2: The degree of tangential bias for dihadron con'elations in the ZOWW formulation of energy loss f?] and in the 
ASW formulation [6J. 



for fluctuations around this average, and as a result the correlation between position and energy 
loss is weakened. More differential observables, such as back-to-back correlations as a func- 
tion of the angle with the reaction plane, which are directly sensitive to the degree of tangential 
bias, have been calculated |,6] and are currently being measured [8]. They may ultimately also 
probe the steepness of the initial medium density profile and allow to make a distinction between 
Glauber and CGC initial conditions in hydrodynamical evolution models and hence also have 
relevance for the discussion of the viscosity of the medium. If this comparison between theory 
and experiment can be made quantitatively, it would be the first real success of tomography in 
the sense of measuring a density distribution. 

2.2. Fully reconstructed jets 

If one is interested in studying correlations on the same side of a hard trigger hadron, one 
probes subleading shower hadrons. Continuing to correlate more and more particles down to 
lower momentum scales, one eventually includes all hadrons of a shower in the analysis and 
hence arrives at fully reconstructed jets. High Pj- physics with jets has a number of conceptual 
advantages over physics with single hard hadrons. First, jets in which the original shower initiator 
momentum is shared across many soft hadrons are frequent, those in which a single hard hadron 
takes most of the initial momentum are not. Thus, jet physics has higher statistics and can reach 
out to larger partonic pr than single hadron physics. Second, the shower evolution in the medium 
happens against a 'meter stick' and a 'clock' given by the medium extension and lifetime, which 
allows to probe the jet evolution not only in momentum space but also in position space. But 
most important, while different models in the energy loss picture may agree on the amount of 
energy lost from the leading parton, disucssing fully reconstructed jets also requires a picture 
how this lost energy is redistributed. In particular, different pictures of leading parton energy 
loss lead to characteristic effects in the distribution of subleading hadrons. This was pointed out 
first in 121 where it was argued that radiated gluons in a radiative energy loss picture lead to a 
charcteristic enhancement in the low f j- region of the longitudinal momentum distribution inside 
medium-modified jets as compared to jets in vacuum. 

Theoretically, the most promising approach to medium-modified jet physics are Monte-Carlo 
(MC) simulations. They build on well-tested vacuum shower evolution codes such as PYSHOW 
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Figure 3: Nucleai' suppression factor Raa computed for a medium-modified sliower in an analytical approximation |^ 
(left) and in Q-PYTHIA El (right) compai'ed with PHENIX data O. 



Ilioll or HERWIG under the assumption that only the partonic evolution of the shower 
is affected by the medium, but that hadronization happens sufficiently far outside the medium 
that it can be treated as in vacuum. Such MC codes have the advantage that they treat energy- 
momentum conservation explicitly at each vertex and furthermore they allow to analyze the out- 
put with the same techniques used to find jets in experimental data. Currently, three major codes 
are being developed: JEWEL follows the philosophy to modify the interaction of shower 

and medium using perturbative QCD under the assumption that the medium constituents can be 
resolved as thermally distributed quarks and gluons. Q-PYTHIA and Q-HERWIG ifTiiflsIl are 
the direct extension of the ASW formaUsm yj from leading parton energy loss to the evolution 



of a whole shower. Finally, YaJEM Hl6l I17L 11811 is based on modifications of parton kinematics 
in the shower where the medium eff'ect is introduced in an ad hoc way in different scenarios (in- 
duced radiation and drag force) to avoid any a priori assumptions about the nature of the medium 
degrees of freedom. Jet observables like the longitudinal distribution of hadron momenta in the 
jet, the thrust distribution or the angular distribution of hadrons in the shower have been obtained 
in all codes. However, the medium-modified shower picture should also account for observables 
sensitive to the leading shower hadron only. 

As an example, Fig.[3]shows the nuclear modification factor Raa computed in the analytical 
approximation of [9] for a schematic treatment of the medium and in Q-PYTHIA averaged over 
a cylindrical medium geometry [14]. Similar results have also been obtained in YaJEM averaged 
over a 3-d hydrodynamical model llldllSH . In all cases, a falling trend of Raa with Pj is observed, 
regardless of the averaging over the medium. At present this is not supported by the trend of 
the data which seems to be increasing, although errors are substantial. However, should the 
increasing trend of the data be confirmed by a higher statistics analysis, it would indicate a 
conceptual issue which is not properly understood. In a similar line, a detailed investigation of 
the LPM eff'ect in the MC [13] indicates that the dependence of radiative energy loss is only 
visible for a certain set of assumptions. If these are relaxed to a more realistic scenario, the 
Lr dependence is hardly visible. However, phenomenologically this poses a problem, as hard 
back-to-back correlations strongly favour a parametric L? dependence of energy loss. It appears 
necessary at this point to test all in-medium shower codes thoroughly against existing leading 
hadron observables before they are applied to predict jet physics. 
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3. The Cone and the Ridge — medium recoil 

While in-medium shower evolution codes treat the production of soft partons perturbatively 
and assume their hadronization can be described as in vacuum, there is no reason to assume that 
this corresponds to the actual dynamics of a jet in the medium at low Pj which is in all likelihood 
non-perturbative. Instead, observed correlations of low Pj hadrons with a hard trigger may rather 
reflect the way the medium responds (in terms of recoil or Shockwaves) to the energy deposition 
by a propagating jet. Thus, more phenomenological models are currently used in attempts to 
describe the soft correlations on the same side and the away side of a trigger hadron. 

In discussing the general structure of the correlation pattern, it is useful to keep a few things 
in mind. If one triggers on a high Pj hadron on the near side, there must be a correlation on the 
away side due to momentum conservation. The interesting physics is only in the dependence of 
this away side correlation on the azimuth 0, the rapidity 77 and Pj, in other words in the precise 
way energy and momentum are balanced. In addition, other physics may lead to a correlation on 
the near side, for example the leading hadron in a jet is correlated with subleading hadrons, but 
this correlation is not required by momentum conservation but arises due to QCD. In a medium, 
other phenomena may also be acidentially correlated with a hard hadron. For example, hard 
hadrons tend to emerge perpendicular to the medium surface because this minimizes the in- 
medium path and hence energy loss. At the same time, radial flow is perpendicular to the surface 
because this is the direction of the pressure gradients. But although this correlates flow with the 
trigger, the two phenomena are unrelated in the sense that radial flow is in no way be caused 
by the hard hadron. Finally, any near side phenomenon correlated with the trigger (be it an 
accidential correlation or not) will also cause a correlation on the away side due to momentum 
conservation. All these effects need to be disentangled carefully. 



3.1. The Ridge 

The most prominent near side correlation associated with a trigger is the so-called ridge, 
a structure which is visible at low Pj, extends out to at least 4 units of rapidity in rf with an 
almost unchanged width in [20.1 . These features practically rule out any connection of the ridge 
correlation with energy loss. The ridge cannot contain radiated quanta, as Raa places a constraint 
on the average amount of energy lost into the medium. A number for the lost energy, given that a 
trigger was observed, has been obtained e.g. in the formalism of [2 1.1 and is (9(0.5) GeV. This is 
insufficient to explain the strong hadron production out to 4 units in rj which takes a substantial 
amount of energy. An alternative explanation of the ridge as a medium recoil due to energy lost 
into the medium fails, since due to the coherence time effect, energy deposition into the medium 
happens comparatively late, but causality requires large spacetime rapidity correlations to arise 
early. However, if one invokes medium recoil, the medium spacetime rapidity must equal its 
momentum rapidity due to the observed Bjorken flow of the medium. Thus, the ridge cannot 
be interpreted as a medium recoil. Besides, a ridge is also observed in untriggered correlations, 
indicating that its correlation with a hard trigger is accidential. On the other hand, the properties 
of the ridge seem well accounted for in models which view the ridge as an initial state such as 
the flux tube model 1 22i |23 ] . 



3.2. The Cone 

The away side correlation at low Pj takes the form of a broad double-hump structure in 
(its structure in 77 is not easily observed as the rapidity of the away side parton is only weekly 
constrained by kinematics). It is most commonly interpreted as conical emission. Early models 
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to explain the cone in terms of Cherenkov radiation 12411 predicted a Pj dependence of the cone 
opening angle which was not found in the data. Therefore the current theory efforts to describe 
the cone focus on the hydrodynamical response to energy deposition into the medium. 

The main conceptual issue is finding the source term for the hydrodynamical evolution which 
describes the physics of energy loss into the medium adequately. Such efforts include a Bethe- 
Bloch source term [25], a source term based on ASW energy loss [26J, or the Higher Twist 
formulation of energy loss ll27ll and also ideas utilizing the AdS/CFT correspondence which 
descibes the physics of a gauge theory in the strong coupling limit |*28i l29ll . 

A problem of the theoretical descriptions of the cone is that many different scenarios man- 
age to result in a double-hump structure of the away side correlation pattern, for example bulk 
statistical momentum conservation and elliptic flow ifsoll . an energy loss source term and hy- 
drodynamical evolution [26], recombination in a linearized hydrodynamical model [31] and the 
recoil of the bulk medium to an initial state flux tube correlation 13211 . This, and the observation 
that the observed angle is not even different at SPS energies makes it likely that the observed 
double-hump structure does not reflect details of the dynamics, but is rather driven by a combi- 
nation of energy-loss caused surface bias, a position-flow correlation in hydrodynamical models 
and a kinematic bias towards alignment of radial flow and large angle correlation, relatively 
independent of the underlying dynamics as long as there is a generic tr ansp ort of energy and mo- 
mentum to large angles. Such a scenario has been explored in lf33i[34 , 35 1 and found qualitative 
agreement with the data on 2-particle correlations. There are also measurements of 3-particle 
correlations, however here neither the experimental analysis nor the theoretical calculation are 
conclusive yet. 



4. Conclusions 

In discussing correlations with respect to a hard hadron trigger as a tomographic tool, it is 
useful to make a distinction between hard and soft associate hadrons. The reason is that due to 
the presence of a hard scale, the virtuality, a hard process in medium and the subsequent partonic 
evolution of a shower remain governed by perturbative QCD dynamics, even if the medium itself 
and its effect on the shower evolution is not. This argument, however, does not apply to soft 
correlations in the typical momentum range of the bulk — here, one cannot maintain an easy 
distinction between recoiling medium and modified shower. 

In the hard sector which can be described in the leading parton energy loss picture, highly 
constrained calculations are possible and being done and begin to answer questions about the 
medium. In practice, this requires to give up schematic descriptions of the medium (Bjorken 
cylinders) and include all constraints of bulk observables by modelling the medium with e.g. 
state of the art hydrodynamical models, and in addition to apply a model not only to single 
hadron supression but also to more differential hard observables, including correlations. Such 
calculations have established reasonably well that elastic energy loss is a small contribution in the 
light quark sector and that the core of the medium is not completely black. On the experimental 
side, more differential high Pj observables are reaching the precision necessary to distinguish 
between models. For observables sensitive to subleading hadrons in the shower, and eventually 
fully reconstruced jets, MC codes for in-medium parton showers are currently being developed 
and tested. They investigate various different ideas how to implement the effect of the medium. 
However, for most codes no thorough comparison with leading hadron observables has been 
done yet, and there are hints of unresolved issues, such as the decreasing trend of Raa with Pj 
or the almost linear growth of energy loss with pathlength. 
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Soft correlations have to be modelled in more phenomenological approaches. There is 
mounting evidence that the near side ridge correlation is not directly connected with energy loss. 
However, on the away side, due to momentum conservation it is quite clear that at least part of 
the cone correlation must have such a connection. At present, the data favour models describing 
the cone as a recoil of the bulk medium over those describing it by a medium modified fragmen- 
tation process. Various ideas for the recoil dynamics are explored, however there are indications 
that the observed correlation is mainly driven by a combination of different biases rather than 
the underlying dynamics. Despite good experimental data, a consistent theory of high Pj bias 
and flow-modified medium recoil is very complicated and at present only capable of qualitative 
answers. Any claims that one would already have measured medium properties such as the speed 
of sound appear hence premature. 
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